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Abstract 

The nuclear j3~ decay from the ground and some excited states of the three-electron 8 Li and 
9 Li atoms is considered. The final state probabilities for the arising Be + ion are determined 
numerically with the use of highly accurate bound state wave functions of the Li atom and Be + 
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ion. The probability of electron ionization during the nuclear /3~ decay of the Li atom is evaluated 
numerically. We also discuss a possibility to observe the double j3~ decay by using the known 
values of the final state probabilities for the regular nuclear f3~ decay. 
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I. INTRODUCTION 

In our earlier study [l| we considered the atomic excitations during the nuclear /3~-decay 
in light atoms and ions. Here we want to study the nuclear /3~-decay in the three-electron 
Li-atom. The main attention below is given to the 8 Li and 9 Li atoms which are of interest 
in some industrial applications. In general, the /3~-decay of the Li atom(s) can be written 
in the following form 

Li^Be+ + e~ + /7 (1) 

where the notation e~ stands for the fast electron emitted during the nuclear j3— decay, 
while V designates the electron anti-neutrino. In general, the nuclear /3~-decay of the Li 
atom leads to the following re-distribution of the bounded atomic electrons. In the result 
of such a redistribution of the incident electron density the final Be + ion can be found in 
a variety of bound states, or even in a number of unbound states. Briefly, this means the 
formation of the two-electron Be 2+ ion during the nuclear (5~ decay of the Li atom. In actual 
applications it is important to predict the probabilities to form the final Be + ion in different 
final states. Note that there are a few selection rules which are applied to the /3 — decays 
in atomic systems (see, e.g., Q). Briefly, these rules can be formulated as the conservation 
laws for the angular momenta L and the total electron spin S. In addition to this, the wave 
functions of the incident and final system must have the same spatial parity. For instance, 
if the incident Li atom was in its 3 2 P(L = 1)— state, then the final Be + can be found only in 
one of its n 2 P(L = l)states. In other words, after the nuclear /3— decay of the Li atom in the 
3 2 P(L = 1)— state it is impossible to detect the final Be + ion, e.g., in the 3 2 S(L = 0)— state, 
or in the 3 2 D(L = 2) -state. 

The advantage of considering three-electron atoms and ions is obvious, since the wave 
functions of such systems can be approximated to very good numerical accuracy. For sim- 
plicity, everywhere below in this study we shall assume that the original Li atom was in its 
ground 2 S(L = 0)— state. The choice of the ground state of the incident Li-atom is not a 
fundamental restriction for our method. Formally, such a state can be arbitrary, e.g., either 
ground state, or 'vibrationally', or 'rotationally' excited atomic state with the given angular 
momentum L (see below) and the total electron spin S. 

By analyzing the properties of the known Li-isotopes one finds that there are two f3~- 
decaying isotopes of lithium: 8 Li (rg ~ 0.84 sec) and 9 Li (rg ~ 0.17 sec). These two 
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isotopes are formed in the (n; 7) — reactions during thermonuclear explosions in which light 
thermonuclear fuel ( 6 LiD) is compressed to very high densities p > 100 g-sm~ 3 by extremely 
intense flux of soft X— ray radiation from the primary. Larger compressions mean, in general, 
the larger output of these two lithium isotopes. In the laboratory, the 8 Li isotope is produced 
with the use of the (n; 7)— reaction at 7 Li. In contrast with this, the 9 Li isotope is produced 
by using either (d;2p) — and (n;p)— reactions with the 9 Be nuclei, or (t; p)— reaction with 
the nuclei of 7 Li. 

The life-times of these two lithium isotopes are relatively short from the chemical point 
of view. Therefore, it is hard to study the regular chemical properties of these isotopes. 
An alternative approach is based on detail analysis of the optical radiation emitted by the 
final Be + ions which are formed after the nuclear f5~ decay of these two Li-isotopes. This 
can be achieved, if we know the correspoding final state probabilities, i.e. the probabil- 
ities to form the final Be-ions in certain bound states. The first goal of this study is to 
evaluate the final state probabilities of formation of various final states in the arising Be + 
ions. Note that all evaluations of the final state probabilities duringthe nuclear /3~-decay in 
atoms and molecules are based on the sudden approximation {3}, Q which applies to both 
atomic systems (original atom and final ion) involved in the process. The sudden approx- 
imation is appropriate for all fi^— decaying atoms, since the velocities of the f3 ± electrons 
are significantly larger than the velocities of regular atomic electrons. 

The final state probabilities, i.e. probabilities to form different atomic states during 
nuclear /3~-decay in various light atoms, have been evaluated numerically in a number of 
earlier papers (see, e.g., [l], Q). All such evaluations, however, have been based on the 
assumption that the total number of bounded electrons does not change during the nuclear 
decay. In reality, the nuclear /3-decay in light atoms often leads to an 'additional' 
electron ionization. For the Li atom this process can be written in the form 

Li Be 2+ + e~ + /T + V (2) 

where e~ designates the secondary atomic electron which becomes free during atomic f3~- 
decay. It is interesting to evaluate the probability of this process and obtain the actual 
energy spectra of the emitted secondary electrons. Formally, all secondary electrons emitted 
during atomic /3-decay must be considered as the 5— electrons. On the other hand, the 
original definition of 5— electrons means that these electrons are fast and their total energies 
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significantly exceed the usual energies of 'regular' atomic electrons. The energy of the 
free electron from reaction, Eq.flSJ), is comparable with atomic energies. Therefore, here 
we deal with the regular atomic ionization during j3~ — decay. In earlier works the process 
of additional ionization only from the atomic K— shell was considered (see discussion and 
references in |5j|). 

The main goal of this study is to determine the final state probabilities to form various 
bound states in the Be + ion. These calculations are discussed in the fourth Section. Another 
aim of our study is to evaluate the probability of 'additional' ionization during the nuclear 
/3~-decay and investigate the energy spectrum of secondary electrons emitted during the 
nuclear /3~-decay. This problem is considered in the third Section. We also briefly investigate 
the old standing problem of the double nuclear /3— decay. Concluding remarks can be found 
in the last Section. 

II. EVALUATION OF THE FINAL STATE PROBABILITIES FOR THE BOUND 
STATES. 

As follows from the general theory of perturbations in Quantum Mechanics (see, e.g., 
j^]) in sudden approximation the final state probabilities are determined as overlap integrals 
between the wave function of the incident atomic system (i.e. the wave function of the Li 
atom in our case) and the wave function of the final atomic system (i.e. the wave function of 
the Be + ion). To compute such an three-electron integral we need to assume that the total 
numbers of bound electrons in the incident and final atomic systems are the same. In sudden 
approximation the general formula for the transition probability Wif for the transitions from 
the incident z-state into the final /—state takes the form (see, e.g., jsj) 



where Vfi is the overlap integral computed with the use of time-independent incident and 
final atomic wave functions, i.e. 



where \&Li(xi, Xg, X3) and ^ Be +(xi, X2, X3) are the wave functions of the Li-atom and Be + 
ion, respectively. The derivation of the formula, Eq.flS]), is based on the facts that: (1) 





Vfi = (^Li(xi,X 2 ,X 3 ) I ^Bc+(X1,X 2 ,X 3 )) 



(4) 
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the velocity of the (3— electron is substantially larger than the velocities of atomic electrons, 
and (2) the final ion does not move during the nuclear ft^— decay. In atomic units we have 
Ti = l,m e = 1 and e = 1 and, therefore, Wif =\ Vfi | 2 . The notation Xj in Eq.(jl]) designates 
the spin-spatial coordinates of the i— th electron, i.e. Xj = (rj,Sj). Note that in some works 
the integral Vfi, Eq.(j3J, (or the ratio ^P-) is called the probability amplitude. The two wave 
functions in Eq.pj depend only upon spatial and spin coordinates of three electrons and do 
not depend upon the time t. All wave functions used in Eq.(Tj]| are assumed to be properly 
symmetrized in respect to all spin-spatial permuations of identical particles (electrons). 

As follows from Eq.flU) the final state probability for the /3~-decay in the 8 Li and 9 Li atoms 
can be determined, if the wave functions of the incident and final atomic systems (bound 
states) are known. The construction of highly accurate variational wave functions for three- 
electron atoms and ions is considered in the fourth Section. The final state probabilities 
determined with the use of such wave functions can be found in Tables I and II. Here we 
assume that the incident Li atom was in its ground 2 S(L = 0)— state (2 2 S— state). It should 
be mentioned that in reality the incident Li atoms are formed in the (n; 7)— , (n; p)— and 
some other nuclear reactions with neutrons of different energies (see above). In such cases it 
is hard to expect that all incident Li atoms will always be in the ground 2 2 S(L = 0)— state. 
In fact, these /^"-decaying Li atoms can be found in a variety of the rotationally and/or 
vibrationally excited states. Very likely, after reactions with neutrons the incident Li atom 
before nuclear /?~-decay will move with the non-zero speed in some direction. Therefore, 
some other (excited) bound states in the Li atom must also be considered as the incident 
atomic states before the nuclear /3~-decay. 

Numerical computation of the overlap integrals, Eq.(jl]), is reduced to calculations of some 
separated integrals, which include different spin components of the incident and final atomic 
wave functions. For instance, let us discuss the construction of three-electron variational 
wave function of the Li atom. Without loss of generality, below we restrict ourselves to the 
consideration of the ground 2 S(L = 0)— state of the Li atom. As is well known (see, e.g., 
jf|, {7 1) the accurate variational wave function of the ground (doublet) 2 S{L = 0)— state of 



the Li atom is written in the following general form 

*(Li) L=0 = if}L=o(A; {rij})(a{3a - f3aa) + (f) L =o(B; {r ij })(2aafi - f3aa - afia) (5) 
where i(jl =0 (A] {r^}) and (j)L =0 (B; {r^}) are the two independent radial parts (= spatial 
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parts) of the total wave function. Everywhere below in this study, we shall assume that all 
mentioned wave functions have unit norms. The notations a and (3 in Eq.flS]) stand for the 
one-electron spin-up and spin-down functions, respectively (see, e.g., |8|). The notations A 
and B in Eq. §5§ mean that the two sets of non-linear parameters associated with the radial 
functions i/) and (b can be optimized independently. In general, each of the radial basis 
functions in Eq.(j5]) explicitly depends upon all six interparticle (or relative) coordinates 
r i2, ^13, ^23, ^14, ?"24, ^4, where the indexes 1, 2, 3 stand for the three electrons, while index 4 
means the nucleus. 

The actual atomic wave function in an atomic system must be completely antisymmetric 
with respect to all electron spin-spatial variables. For three-electron wave function this 
requirement is written in the form ^4 123 \I / (1, 2, 3) = —^(1,2,3), where the wave function 
\l/ is given by Eq.(j5]) and A e is the three-particle (or three-electron) antisymmetrizer A e = 
e— P\2— P13— A3+A23+ A32- Here e is the identity permutation, while Aj is the permutation 
of the 2-th and j-th particles. Analogously, the operator P^ is the permutation of the z-th, 
j-th and /c-th particles. 

Suppose that the incident three-electron wave function of the Li atom has been chosen 
in the form of Eq.([5]). By applying the antisymmetrizer ^123 to the first part of the total 
wave function, Eq.(j5]), one finds 

-4.123 4>l=o(A; {rij})(a(3a — (3aa)} = (e?/>) (ct/fo — /3aa) + (Pi2ip)(a^a — fiaot) 
— {Pi?,ip){ctf3a — aa/3) — (p23ijj)(aa(3 — (3aa) + (A23VO — aj3a) 

+ (-Pi32^)(/3aa - aa/3) (6) 

where the notations (Piji/j) and (Ajfc'*/') mean the permutation operators which act on the 
coordinate functions only. Analogously, for the second part of the total wave function one 
finds 

vAi23 <Pl=o(B'i {rij})(2aa/3 — j3aa — a/3a)] = {ecj)){2aaj3 — (3aa — af3a) 
— ( (2aa/3 — (3aa — a/3a) — (Pi3<j))(2f3aa — aa[3 — aj3a) 
— (p23(f>)(2a(3a — aa/3 — f3aa) + (Pi23(b)(2(3aa — afia — aaf3) 

+(A320) (2a/3a — aa/3 — j3aa) (7) 

where the notations (Aj0) an d (Ajfc</>) mean the permutations of the spatial coordinates in 
the (j)L=o(B] {ry}) radial function, Eq.([5]). 



Now, by using the expressions, Eqs.flH]) and ([7]), we can obtain the formulas which can 
be used in computations of the final state probabilities in the case of the nuclear /3~-decay, 
Eq.([T]), in the three-electron Li atom. For instance, if the final wave function has the same 
spin-symmetry, i.e. it is written in the form 

tyfi = ipfi(ri, r 2 , r 3 )(a/3a - (3aa) + <^(ri,r 2 , r 3 )(2aa/3 - (3aa - a /3a) (8) 

then the final state probabilities are determined with the use of the following formulas 

P^ = (^(n, r 2 , r 3 ) | ^(2e + 2P 12 - P 13 - As - P 123 - P lm )^u{A- {n 3 })) (9) 

P H = (^(n, r 2 , r 3 ) | ~(Aa - As + A23 - Aaa)^^; M)) (10) 
P^ = (^(ri.^.ra) | ~(A 3 - As + A23 - A32WB; {r^})) (11) 



P 00 = (^-(r^ra.ra) | — =(2e - 2P 12 + P 13 + F 23 - F 123 - P 132 )<t>u(B; {r^})> (12) 

a 

Note that these formulas coincides with the known formulas [9J which correspond to the 
doublet — > doublet transition in the three-electron atomic systems. Briefly, this means 
that the both incident and final atomic states contain three-electrons in the doublet spin 
configuration (the total electron spin equals |). This case corresponds to the 'classical' 
/3 ± -decay in few-electron atoms, when the incident and final electron configurations has the 
same L and S quantum numbers. 

In reality, another process is also possible in few-electron atoms during the nuclear (3 ± - 
decay in few- and many-electron atoms. This process leads to the formation of the final 
ion/atom in some excited states. For instance, consider the case when the three final elec- 
trons form the doublet configuration with the spin function aa/3. It is clear that such 
a wave function cannot represent the ground state of the Be + ion. However, some ex- 
cited states (with vacancies in the internal electron shells) can have this spin function, e.g., 
Is2s3p— , Is2s4d— and Is3s3p-states of the Be + ion. Another example is discussed in the 
next Section. It represents an additional electron ionization during the nuclear j3~ decay in 
three-electron atom. If this free electron moves away in the /3— spin state, then the final Be 2+ 
ion can be only in its triplet spin state (not singlet state). Formally this means formation 
of the final ion in an excited state (with some vacancies on its internal electron shells). In 
this case the formulas for the final state probabilities take the form 

Ptm> = (^(ri,r 2 ,r 3 ) | ^(As - As + A23 - As 2 )^Li(^ {r t ,})) (13) 
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P t r<t, = (^(ri, r 2 , r 3 ) | ^=(2e - 2P 12 + As + As - A23 - Pin) MB; {r,,})) (14) 

where it is assumed that the incident electron wave function was written in the form of 
Eq.flH]). These formulas indicate clearly that the probabilities to find the final Be 2+ ion in 
the excited triplet spin states are not zero. In all earlier studies the transitions to the final 
atomic states with different spin states were ignored. Moreover, any possibility to form the 
final few-electron ion/atom in excited states (with some vacancies on internal electron shells) 
during the nuclear /3 ± -decay was rejected. The actual existence of such transitions is a great 
achievement of this study. 

In general, during the nuclear /3~-decay of the Li atom the final Be + ion can be formed in 
many different bound and/or unbound states. If such a state is unbound, then we deal with 
the additional ionization during atomic j3~ decay. It is discussed in the next Section. This 
process is of great interest, since it often leads to the formation of the final ion in an excited 
state(s) with various vacancies in the internal electron shells. For light atoms and ions this 
means a possibility to observe emission of the optical quanta after the nuclear (3^— decay in 
many-electron atoms with the total number of electrons > 3. 

III. ADDITIONAL IONIZATION DURING (3 -DECAY. 

The probability of ionization (or 'additional' ionization) of the final Be + ion during the 
nuclear j3— decay can also be evaluated with the use of the sudden approximation. In this 
case the final wave function is constructed as the product of the bound state wave function of 
the two-electron Be 2+ ion and the wave function of the unbound electron which moves in the 
central Coulomb field of this two-electron ion. To determine the corresponding final state 
probability one needs to compute the following overlap integral between the wave functions 
of the incident and final atomic systems 

Afi = (^Li(xi,X 2 ,X 3 ) I ^Bc+( X 1 1 X 2)0(X 3 )) (15) 

where 0(xs) is the wave function of the unbound electron which moves in the Coulomb field of 
the Be 2+ ion. This function must include the continuous parameter k which is the electron's 
wave number (see below). It should be mentioned that such an 'additional' ionization has 
nothing to do with the interaction between the emitted (3~ electron and atomic electron. In 
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fact, the additional ionization is directly related with the presence of the non-zero component 
~ 0(xa) in the incident atomic wave function. 

The probability of additional ionization has been determined for a number of 

n n 

(3 —decaying atoms in a number of earlier studies (see, e.g., [4| and [5[). These works, how- 
ever, were restricted to the analysis of the electron ionization from the internal K— shells 
only. In this case the original problem was reduced to the solution of the model one-electron 
problem. Analogous reduction for few-electron atomic systems is much more difficult to per- 
form, since all electron-nucleus and electron-electron coordinates are not truly independent. 
It complicates accurate computation of integrals which involve such coordinates. Neverthe- 
less, numerical computations of the final state probabilities can be conducted even with the 
use of highly accurate wave functions known for many few-electron atoms. In this Section 
we discuss some details of such calculations. 

In atomic units the explicit form of the one-electron wave function is </>(r) = <pki(r)Yi m (n) , 
where 4>ki{f) is the one-electron radial function, while YJ m (n) is the corresponding spherical 
harmonics and n = ^ is the unit vector accosiated with r. In this Section the parameter k 



is k = y ™$ E = V2E (in atomic units). The explicit formula for the radial function 4>ki{r) 
(in atomic units) is (see, e.g., Q) 

C 

Mr) = j^Yy (2Qkr) 1 ■ ex P (-tQkr) ■ 1 F 1 (— + 1 + 1,21 + 2,2iQkr) (16) 



where iFi(a, b;x) is the confluent hypergeometric function (see, e.g., [U]), while Cm is the 
following constant 



C 



kl 



InQk 



i-«p(-&)j s x = y\r 1 q 2 p 



IW* 2 + 7^ ( 17 ) 



Qk> s=l 

In these two equations the parameter Q is the electric charge of the remaining double-charged 
(positive) ion, i.e. Q = 2. In reality, this parameter must slightly be varied (around 2) to 
obtain better agreement with the experimental data. Such variations formally represent 
ionizations from different electronic shells of the incident Li atom. 

Accurate numerical computations of the final state probabilities during the nuclear /3~ 
decay in few-electron atoms with additional electron ionization are very difficult to perform, 
since all highly accurate wave functions of the bound states explicitly include the electron- 
electron coordinates. As a rule, the better accuracy of the bound state wave function 
means more complete and accurate involvement of the terms which describe various electron- 
electron correlations. On the other hand, the crucial step of the whole procedure is the 



numerical and/or analytical computations of the Fourier transformation of the one-electron 
wave function. This corresponds to the free motion of the final electron. During that 
step of the procedure (Fourier transform) it is better to consider all electrons as particles 
independent of each other, i.e. ignore all electron-electron correlations. In the general case, 
this two-fold problem has no simple solution which is accurate and relatively simple for 
Fourier transform at the same time. 

In this study we have developed an approximate procedure which can be used to perform 
approximate numerical evaluations for the /^"-decaying isotopes of the three-electron atoms. 
In this approach the trial wave function is constructed as the sum of many terms and 
each of these terms contains the products of the electron-nucleus functions. None of the 
three electron-electron coordinates r 32 ,r 31 ,r 2 i is included in such trial wave functions. For 
the ground state (the doublet 2 S(L = 0)— state) of the Li atom the radial wave function 
^L=o(A; {rij}) is chosen in the following form: 

ipL=o{ru, r 24 , r 34 , 0, 0,0) = ^ C k r^ {k) r^ k) r^i {k) exp(-a k r 1A - (3 k r 2 4 - 7fc r 3/t) (18) 

fe=i 

= £ C , fc rJ niW rJ l2( *Vj ,a( * ) exp(-afcri - (3 k r 2 - 7 fe r 3 ) 

k=l 

where C k are the linear (or variational) coefficients, while mi(k),m2{k) and m 3 (/c) are the 
three integer (non-negative) parameters, which are, in fact, the powers of the three electron- 
nucleus coordinates = (i = 1, 2, 3). Below, we shall assume that the trial wave 
function Eq.( fl8l) has a unit norm. Furthermore, in all calculations performed for this study 
only one spin function Xi{Xi = a P a ~ Pcua) is used. It is clear that the wave function 
Eq. ({TBI) contains only electron- nuclear coordinates and does not include any of the electron- 
electron coordinates. The real (and non-negative) parameters a k ,P k ,'j k are the 3N S varied 
parameters of the variational expansion, Eq. fllSI) . The wave function, Eq. fllSp . must be 
properly symmetrized upon all three electron coordinates. 

The principal question for the wave function, Eq. (1181) . is related to its overall accuracy. 
If (and only if) such an accuracy is relatively high, then the trial wave function, Eq. (fl8l) . 
can be used in actual computations of the probability amplitudes. For this study we have 



constructed the 23-term variational wave function shown in Table I of Ref . jlOj . This wave 
function is represented in the form of Eq. ffTBT) and contains no electron-electron coordinates. 
All sixty nine (69 = 3 x 23) non-linear parameters in this trial wave function have been 
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optimized carefully in a series of bound state computations performed for the ground 2 S(L = 
0)— state of the Li atom. Finally, the total energy E of the ground 2 S— state of the °°Li 
atom obtained with this independent-electron wave function is -7.44859276608 a.u. Note 
that such an energy E is close to the exact total energy of the ground state of the °°Li atom. 
It indicates a good overall quality of our approximate wave function with 23 terms which 
does not include any of the electron-electron coordinates r^^r^^r^. This wave function is 
used in computations of the final state probabilities (see below) for the nuclear (3~ decay 
with additional electron ionization in the three-electron Li atom. 

Note also that in atomic physics based on the Hatree-Fock and even hydrogenic ap- 
proximations the ground state in the Li atom is designated as the 2 2 S— state, while in the 
classification scheme developed in highly accurate computations the same state is often des- 
ignated as the l 2 ,?— state. This classification scheme is very convenient to work with trully 
correlated few-electron wave functions. It is clear that no hydrogenic quantum numbers are 
good in such cases, and we have to use the more appropriate (and convenient) classification 
scheme. To avoid confllicts between these two classification schemes in this study we follow 
the system of notation used earlier by Larsson [7] which designated this state in the Li atom 
as the 'ground 2 6'-state'. 

The wave function of the final two-electron Be 2+ ion arising during the nuclear /3~-decay 
with the additional ionization can also be approximated with the use of basis functions which 
depend upon the electron-nuclear coordinates only and do not include the electron-electron 
coordinate r 2 i. For the bound S(L = 0)— states of the Be 2+ ion such an expansion takes the 
form 

^L=o(n, r 2 , 0) = C k r™ l{k) r™ 2{k) exp(-a k r 1 - j3 k r 2 ) (19) 

k=l 

The use of the approximate wave functions Eqs. tjTBl - ( TL91) which do not depend upon 
explicitly the electron-electron coordinates simplify numerical computations of all integrals 
required for the numerical evaluation of the final state probabilities during the nuclear /3~- 
decay in the three-electron atoms and ions. The remaining part of the problem is the 
analytical computation of the integral between the product of the factor r m exp(— jr) and 
radial function from Eq.f fTBl . Such an integral is computed with the use of the formula (see, 
e.g., Eq.(7.522.8) from (llj): 

r +co T(u + 1) c 

/ exp(-Xx)x u 1 F 1 (a,b;cx)dx = -± r-zF^a, v + 1; 6; -) (20) 

Jo X u + 1 A 
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where T(x) is the usual 7— function (see, e.g., Section 8.31 in [ll|). Our results of numer- 
ical computations of the final states probabilities for the nuclear /3~-decay with additional 
electron ionization will be published elsewhere. 



IV. BOUND STATE WAVE FUNCTIONS OF THE THREE-ELECTRON ATOMS 
AND IONS 

As is shown above to determine the final state probabilities during the nuclear /3~-decay 
one needs to know the accurate wave functions of the incident and final atoms and ions. 
In sudden approximation the angular momentum L, electron spin S and spatial parity n 
of the atomic wave function \1/ are conserved during the nuclear /3~-decay. Therefore, all 
approximate wave functions must be constructed as the eigenfunctions of the operators of 
angular momentum L 2 and total electron spin S 2 . In this study we use the variational wave 
functions constructed with the help of the method of Hylleraas-Configuration Interaction 
(Hy-CI). In general, the wave functions of Hylleraas-type expansion rapidly converge to the 



exact wave functions. The Hylleraas-Configuration Interaction wave function 12|, |l3| is a 
linear combination of symmetry adapted configurations $ p : 

JV 

#Hy-CI = E C P%, % = OCL 2 )MpX (21) 
p=l 

where the spatial part of the basis functions are Hartree products of Slater orbitals containing 
up to one interelectronic distance per configuration: 

n 

= r ij IT M r k, 9k, <Pk)- (22) 

k=l 

If v — 1, then the wave function, Eq. (|2ip . corresponds to Hy-CI. In the case when v = 0, it is 
the usual Configuration Interaction (CI) wave function. The higher powers of the electron- 
electron distances r\- can effectively be reduced to the term (or u = 1). Indeed, all higher 
terms v > 1 can be expressed as a product of r^, a polynomial of ri,rj and some angular 
functions. Also, in Eq. (|2ip iV is the number of configurations used in computations. The 
coefficients C v are determined variationally. The operator 0(L 2 ) in Eq. fl21l) projects over 
the appropriate space, so that every configuration is an eigenfunction of the square of the 
angular momentum operator L 2 . A is the antisymmetrization operator upon all electron 
spin-spatial coordinates and \ is the electron-spin eigenfunction. For the lithium atom and 
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three-electron ions one can choose the total spin function in the one-component form, i.e. 
X = (a/3 — (3a)a. 

The basis functions 0& in this work are the s-, p-, d-, and /-Slater orbitals. Since the 
convergence of Hy-CI wave functions is usually very fast, there is no need to use orbitals 
with higher angular momentum. The unnormalized complex Slater orbitals are defined as: 

^(r) = r^ 1 e-°T(fl, V ). (23) 

where the parameter a is the adjustable variable (for each orbital) and Y l m (9, ip) are the 
complex spherical harmonics. The basis sets employed in this work are n — 4, 5, 6 and 7, 
where the basis n — 4 (for short n4) means the orbital set [ls2s3s4s2p3pAp3dAdAf]). With 
all these orbitals from our basis set we have constructed the most important configurations 
of the S(L = 0,M = 0)-, P(L = 1,M = 0)-, and D(L = 2, M = 0)-symmetries. All 



details of construction of the symmetry adapted configurations <3> p of Eq. ( 12TT) can be found 
in Ref.fl. 

The orbital exponents have been optimized for each atomic state of the Li atom and Be + 
ion. A set of two exponents have been used, one for the K-shell and the other for the odd- 
electron of the L-shell. It was kept equal for all configurations. This technique accelerates 
numerical computations. The obtained results are sufficiently accurate for the purposes of 
our investigation. It is clear that to obtain highly accurate energies one needs to apply 
more flexibility in the exponents. It was shown in recent calculations on the lithium atom 



17]. The virial factor x = — -S3- is used to check the quality of the wave function and 



<T) 

guide the numerical optimization of the exponents in the trial wave functions. 

As for a given basis set the number of possible configurations would be too large, we have 
selected the 'most important' configurations according to their contribution to the total 
energy. In our case the selection criteria is an energy contribution > 10~ 8 a.u. with respect 
to the previous configuration. For that, blocks containing all possible Hy-CI configurations 
of the same type have been filtered and the configurations with less energetic weight have 
been thrown out. More details can be found in [l^]. Note that the length of the wave 
functions varies then for every state, and the selected configurations are for every state 
different. As a result, higher excited states must be approximated with the use of the longer 
trial functions. 

For our calculations in this study we have written a three-electron Hy-CI computer pro- 
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gram in Fortran 90 language. The calculations have been done with the use of the quadruple 
precision. The program has been thoroughly checked by comparing results of our numeri- 
cal calculations with the analogous results obtained earlier by Sims and Hagstrom and 



King 16] for the lithium atom. Note that in such calculations we have observed a complete 
agreement. The energy calculations of the ground and excited S-states of the both Li atom 
and Be + ion are shown in Table I, together with their convergency with respect to the basis 
set used. For the ground states of the Li atom and Be + ion an accuracy of (1.4 - 1.9)-10 -6 
a.u. has been achieved with the techniques described in this paper. For the first two excited 
states of both Li atom and Be + ion the accuracy is of about 5-10 -6 a.u.. In the third and 
higher excited states within a given symmetry of the S-,P- and D-states the accuracy is of 
the order of ±1 • 10~ 3 a.u.. For the accurate calculation of these higher excited states it 
would be necessary to introduce different sets of orbital exponents and to increase the orbital 
basis. However, numerical calculations on higher excited states are sparse in the literature. 

For numerical calculations of the amplitudes and transition probabilities during nuclear 
/3~-decays in three-electron atomic systems we have developed the new computer program 
which calculates the overlap integrals, Eq.(|l]), between the wave functions of Li atom and 
Be + ion. The previous step is the calculation of the wave functions of the different states 
of the Li and Be + atoms using the Hy-CI method. The algorithms we have employed for 



the calculation of the kinetic and potential energy integrals can be found in Refs. 
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201 ]. The energy values obtained for the ground and S-, P- and D-excited states are given 
in Ref. fl4|. as so as more details on the calculation and the comparison with the best 
data of the bibliography. For contrary in the calculation of the overlaps we need only 
the overlap-integrals between the configurations, and the coefficients of the Hy-CI wave 
functions. Therefore in this program only integrals of the types (r™ 2 ), (r^r^) are needed, 
while the fully-linked three-electron integral ( ri r 2 2 r 3 13 ) is n °t needed, when the overlap between 
the wave functions containing the terms is calculated. 

In this work we have improved our earlier method of calculations of the final state prob- 

n 

abilities during the nuclear f3 -decay [1]. Now, we calculate the overlap between the wave 
functions of different lenghts. This overlap is the sum of the matrix elements of a rectan- 
gular overlap matrix. This method of calculation has an advantage, since there are several 
possible checks. First, the permutation symmetry of the overlap matrix (^i] 1 !^) = ($2 | ^i) 
and its unit-norm condition, i.e. i^i \ ^>i) = (^2 I ^2) — 1- 
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The convergence of the probability amplitudes and probabilities increases with the im- 
provement of the total energies of the incident and final atomic systems. The final transition 
probabilities are calculated with an error < 0.01% (they are summarized in Table II). We 
have obtained the transition probability from the ground 2 S— state of the Li atom to the 
ground state of the Be + ion 57.71% . The transition probability for the transition from 
the ground 2 S— state of the Li atom to the first excited 2 S— state of Be + ion is ~ 26.51 % and 
to the second excited 2 S— state such a probability is 0.55% . The sum of the probabilities 
is then ps 84.99% and the probability of ionization calculated as 1.0 minus this sum is pa 
15.00% . 

In addition we have computed the transition probabilities from the lower-lying excited 
states of S-, P-, and D- symmetry of the Li atom to the states of the same symmetry states 
in the Be + ion. The results can be found in Tables III, IV and VI, respectively. It is clear 
that the sum of the probabilities of transition from on state of the incident atom to the 
states of the final one must always be less (or equal) unity. In this work we have checked 
this condition everywhere. For instance, the transition probability 2 2 S (Li) — > 3 2 S (Be + ) 
is ps 56.93 %, while the probability of transition 3 2 S (Li) ->■ 4 2 S (Be + ) is « 62.46% . It 
is interesting to note that the probabilities of transition from the incident atom in excited 
and rotationaly excited states are relatively large, showing the possibility of these processes, 
for example the D-transitions l 2 L>(Li) ->■ l 2 D(Be + ) « 37.73% and l 2 L>(Li) ->■ 2 2 D(Be + ) ps 
44.69%. This indicates, clear that from each state in the f3~ — decaying atom there only are 
a few final (partner) states. The probability of ionization (or more properly that the final 
electronic state is a Q 1 S, Q 1 P or Q 1 D or higher) is also calculated. 

Finally, from the methodogical point of view, it is interesting to compare the results 
obtained by using Hy-type wave functions and CI wave functions. Traditionally, in the 
literature on the nuclear /3~-decay calculations in atoms it has been generally assumed 
that Hy-type wave functions are needed, but a real proof of this cannot be found anywhere. 
Therefore in this work we want to prove that fact. In Table VI we show the same calculations 
obtained with accurate CI L-S wave functions using Slater orbitals. The probabilities have a 
poor convergency compared with the ones obtained by Hy-type wave functions. This means 
when improving the accuracy of the energy the calculated amplitude does not converge. 
We have investigated this problem by using the double and quadruple precision, showing 
that there is no loss of precision in the coefficients of the wave functions, but the loss of 
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precision is in the trial wave functions themselves. The reason is that the /?-decay is a nuclear 
phenomena where the change of charge in the nucleus plays an importan role. The different 
performance of Hy-CI and CI wave functions is due to their methodological differences in 
the description of the electronic density near of the nucleus. 

In recent calculations of the transition probabilities for the nuclear /?~-decay of the helium 
atom into the Li + ion (He — > Li + ), Wauters and Vaeck {2l| by using the basis set of B-splines 
obtained the total energies of the incident He atom and final Li + ion which were very close to 
the energies obtained by the Hylleraas wave functions. This is due to the fact that this basis 
set are designed to describe the real shape of the orbitals. However, this is not a standard 
CI method. In general, the transition probabilities during the nuclear /3 ± -decay calculated 
with the CI method are not very accurate. Higher accurate wave functions which describe 
correctly the nuclear-electron cusp and therefore the electron density near of the nucleus are 
needed to calculate precise transition probabilities during various nuclear processes. 

V. ON THE DOUBLE p DECAY. 

The idea of the double nuclear (3 ± — decay in some nuclei was proposed in the middle 
of 1950's by A.B. Migdal and others. Such an assumption was based on observation of 
the nuclear products formed during large-scale thermonuclear explosions. Extremely high 
compressions of the light thermonuclear fuel (mainly 6 LiD) and very intense neutron fluxes 
generated in these explosions, produce a significant amount of neutron-excessive isotopes of 
the heavy transuranium elements. The overall intensity of (3— decays observed in a few first 
seconds after large thermonuclear explosions is very high. However, already after 10 - 15 
seconds the remains of thermonuclear explosive device contain only regular 'slowly '-decaying 
/^"-isotopes of usual heavy elements. The following intensity of the nuclear /3~-decay from 
the products of thermonuclear explosion decreases with time as ~ t -137 [22|. Approximate 
balance between the total number of nuclei of /^"-decaying isotopes which are formed and 
destroyed during large thermonuclear explosions leads to a conclusion that the double 
decay is possible. Originally, the double j3~— decay of neutron-excessive nuclei was proposed 
to explain very high intensities of the nuclear (3— transformations in a first moments after 
large thermonuclear explosions. Formally, there is no restriction (or selection rule) which 
can prohibit the double /3~-decay in atomic nuclei. However, since the middle of 1950's all 
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experimental attempts to detect the double nuclear decay have failed. 

In this study we do not want to discuss neither theoretical significance of the double 
/3~-decay for nuclear physics (or for the finite Fermi systems), nor its possible applications. 
Instead, let us consider the difference in the final atomic probabilities which can be obtained 
in the two following cases: (a) the double nuclear (3— decay, and (b) two consecutive (single) 
nuclear f3~— decays. From atomic point of view we need to compare the time T2/3 for which 
two [5~ particles leave the nucleus with the regular atomic time r e = — n ^ Q2 , where Q is 
the electric charge of the nucleus expressed in e, i.e. Q = Qe. The condition r 2( g <C r e 
means sudden emission of the two fast (5~ particles. In this case the probability amplitude 
is determined as the ovelap integral of the incident and final (atomic) wave functions, Eq.(j3J). 
If the equation of the double (3~ decay is written in the form X — > Z 2+ + /3f + j3 2 + v\ + V2, 
then for the final state probability one finds (in atomic units) 

w fi =| (* x (xi,x 2 , . . . ,x„) I ^z2+(xi,x 2 , . . . ,x n )) I 2 (24) 

In the opposite case, i.e., when r 2/ 3 3> r e , we deal with the two consequent nuclear j3~ 
decays. In this case we need to use the sudden approximation twice. The corresponding 
consequtive equations are X — >■ Y + + /3-f + V x and Y + — > Z 2+ + (3^ +V 2 . The amplitude of 
the final state probability takes the form 

Aft = (^ X (xi,X 2 , . . . ,X„) I ^ Y +(xi,X 2 , . . . ,X n ))(^ Y +(xi,X 2 , ...,X„) | ^Z2+(Xi,X 2 , . . . ,X n )) 

(25) 

and the final state probability is Wfi =\ Afi \ 2 . The sum over all states of the Y + ion will 
lead us back to Eq. (1241) . However, if r 2 ^ 3> r e , then the final state of the Y + ion is uniformly 
defined and Eq.( l25|) can be used in this case only for this unique state of the Y + ion. The 
sum over all 'intermediate' states of the Y + ion is reduced to the one term only. From here 
one easily finds that 

w fi (X -> Z 2+ ) > w fi (X -)• Y + )^(Y+ -)• Z 2+ ) (26) 

In general, the study of the double nuclear /3~-decay in atoms and molecules can be used 
as a natural tool to study electron-nucleus and electron-electron correlations at the femto- 
and attosecond time-scale. Unfortunately, at this moment nobody performs similar research 
even for atoms and ions. 
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VI. CONCLUSION 



We have considered the nuclear /3~-decays in the three-electron 8 Li and 9 Li atoms. The 
final state probabilities to form different bound states in the Be + ion have been determined 
to very good accuracy. 

The Hylleraas-CI wave functions constructed for atoms/ions involved in the /3~-decay 
are substantially more accurate than CI wave functions used in earlier studies. The reason 
is they provide a better description of the electron density near of the nucleus. For the 
first time, the wave functions of the excited states are determined to the same numerical 
accuracy than the wave functions of the ground states. By using such wave functions we 
could determine the final state probabilities to very high accuracy. 

We also discuss a possibility to observe the double nuclear /3~— decay and nuclear 
j3~— decay with the additional ionization. It is shown that the Be 2+ ion formed during 
the last process can be detected not only in the singlet bound states, but also in the triplet 
bound states. It was never observed/predicted in earlier studies on the f3~— decays in atoms. 
Briefly, we can say that the nuclear j3~— decay with the additional ionization in few-electron 
atoms may lead to the fundamental re-structuring of the internal electron shells of the inci- 
dent atom. 
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TABLE I: Convergency of the Hy-CI calculations of the transition probabilities of j3 -decay of the 
lithium atom with respect the basis set employed. Energy in a.u. 



State 


Basis" 


N 6 


Amplitude 


Probability 


Virial 


Energy 


Ref. Ener. c 


Diff. in /ih 


Li 1 2 S 


ni 


308 






2.000004 


-7.4780 5322 25 




7.1 


Li 1 2 S 


n5 


517 






2.000001 


-7.4780 5782 50 




2.5 


Li 1 2 S 


n(> 


620 






2.000000 


-7.4780 5889 25 


-7.4780 6032 39 


1.4 


Be+ 1 2 S 


ni 


308 


0.7596 8128 09 


0.5771 157 


2.000002 


-14.3247 5737 75 






Be+ 1 2 S 


n5 


612 


0.7596 8642 44 


0.5771 234 


2.000001 


-14.3247 6041 19 




2.8 


Be+ 1 2 S 


n(> 


707 


0.7596 8900 66 


0.5771 273 


2.000000 


-14.3247 6139 60 


-14.3247 6317 68 


1.8 


Be+ 2 2 S 


ni 


307 


0.5149 4787 82 


0.2651 713 


2.000000 


-13.9227 5998 01 






Be+ 2 2 S 


n5 


459 


0.5148 9299 62 


0.2651 148 


2.000009 


-13.9227 8162 47 




7.0 


Be+ 2 2 S 


n(> 


535 


0.5149 3042 17 


0.2651 533 


2.000008 


-13.9227 8531 43 


-13.9227 8926 85 


3.9 


Be+ 3 2 S 


ni 


252 


0.0752 1442 67 


0.0056 572 


2.001593 


-13.7985 2045 32 






Be+ 3 2 S 


n5 


372 


0.0742 1482 48 


0.0055 078 


2.000484 


-13.7987 0472 20 




16.3 


Be+ 3 2 S 


n(> 


451 


0.0738 5606 29 


0.0054 547 


2.000154 


-13.7987 0889 23 


-13.7987 1657 


7.7 


Be+ 4 2 S 


n5 


502 


0.0398 3939 53 


0.0015 871 


2.004399 


-13.7347 3172 02 






Be+ 4 2 S 


116 


516 


0.0397 6930 41 


0.0015 816 


2.004126 


-13.7424 8645 38 


-13.7446 306 


3106.7 


Be+ 5 2 S 


n6 


603 


0.0256 1557 46 


0.000 6562 


2.003919 


-13.7152 8280 61 


-13.7162 8624 


1003.4 



"Basis set, i.e. n4 stays for [4s3p2dlf] or [Is2s3s4s2p3p4p3d4d4/]. 
b N is the number of Hy-CI symmetry adapted configurations. 

c Details on the calculations and comparison with reference energies are given in Ref. 14j . 
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TABLE II: Transition probabilities for the nuclear /?~-decay from the ground l 2 «S-state of the Li 
atom into the ground and various excited states of the Be + ion. The probability of ionization from 
Be + ion to Be 2+ ion is calculated as Pi on = 1 — Pi- 



State of Be+ 


Amplitude 


Probability (Pi) 


Pi in % 


1 2 S 


0.7596 8900 66 


0.5771 273 


57.71 


2 2 S 


0.5149 3042 17 


0.2651 533 


26.51 


3 2 S 


0.0738 5606 29 


0.0054 547 


0.54 


4 2 S 


0.0398 3939 53 


0.0015 871 


0.15 


5 2 S 


0.0256 1557 46 


0.0006 562 


0.06 


Total 




0.8499 708 


84.99 


Pion 




0.1500 292 


15.00 
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TABLE III: Transition probabilities for the nuclear (3 decay from the excited 2 2 S-, 3 2 S-, 4?S-, 
5 2 5-states of the Li atom a into the ground and various excited states of the Be + ion. 



States Li — > 


Be+ Amplitude 


Probability 


Pi in % 


2 2 S -> 


1 2 S 


0.2399 9904 24 


0.0575 995 


5.76 


2 2 S -> 


2 2 S 


0.4662 0000 32 


0.2173 424 


21.73 


2 2 S -» 


3 2 S 


0.7575 2492 75 


0.5738 440 


57.38 


2 2 S -> 


4 2 S 


0.0165 6566 49 


0.0002 744 


0.03 


2 2 S -» 


5 2 S 


0.0132 6988 87 


0.0001 761 


0.02 


Total 






0.8492 364 


84.20 


p. 

1 ion 






0.1507 636 


15.08 


3 2 S -> 


1 2 S 


0.1322 9548 78 


0.0175 021 


1.75 


3 2 S -> 


2 2 S 


0.2354 2179 57 


0.0554 234 


5.54 


3 2 S -> 


3 2 S 


0.1207 0198 23 


0.0145 690 


1.46 


3 2 S -» 


4 2 S 


0.8274 1166 41 


0.6846 101 


68.46 


3 2 S -> 


5 2 S 


0.2409 8079 97 


0.0580 717 


5.81 


Total 






0.8301 763 


83.02 


p. 

J ion 






0.1698 237 


16.98 


4 2 S -> 


1 2 S 


0.1024 2628 65 


0.0104 911 


1.05 


4 2 S -> 


2 2 S 


0.1700 6598 04 


0.0289 224 


2.89 


4 2 S -> 


3 2 S 


0.1068 2286 82 


0.0114 111 


1.14 


4 2 S^ 


4 2 S 


0.2356 9455 04 


0.0313 772 


3.14 


4 2 S -» 


5 2 S 


0.0132 6988 87 


0.0001 761 


0.02 


Total 






0.0823 779 


22.52 


p. 

J ion 






0.9176 221 


91.76 


5 2 S -> 


1 2 S 


0.1304 3306 40 


0.0170 127 


1.70 


5 2 S -> 


2 2 S 


0.2194 2766 22 


0.0481 485 


4.81 


5 2 S -> 


3 2 S 


0.1985 1409 22 


0.0394 078 


3.94 


5 2 S -> 


4 2 S 


0.0788 4729 63 


0.0006 217 


0.06 


5 2 S -> 


5 2 S 


0.3681 2183 36 


0.1355 137 


13.55 


Total 






0.2407 044 


24.07 


p. 

J ion 






0.7592 956 


75.93 



°The calculated energy of the first excited S-state of lithium atom 2 2 S is -7.3540 9081 85 a.u. and of the 

higher excited states 3 2 5 -7.3184 8295 19 a.u., 4 2 S -7.3018 6513 00 a.u. an 5 2 S -7.2859 3924 50 a.u. 
6 The calculated energies of the S-states of Be + are the best values of Table I. 
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TABLE IV: Transition probabilities between states of P-symmetry for the nuclear j3 -decay of the 
Li a atom to the Be + ion b . 



States Li — > 


Be + Amplitude 


Probability 


P l in % 


1 2 P - 


1 2 P 


0.6943 3672 59 


0.4821 035 


48.21 


1 2 P - 


2 2 P 


0.6833 4887 87 


0.4669 657 


46.70 


1 2 P - 


3 2 P 


0.0105 6023 91 


0.0001 115 


0.01 


1 2 P - 


4 2 P 


0.0393 4260 79 


0.0015 478 


0.15 


1 2 P - 


5 2 P 


0.0116 1081 46 


0.0001 348 


0.01 


Total 






0.9508 633 


95.09 


p. 

1 ton 






0.0491 367 


4.91 


2 2 P - 


1 2 P 


0.2744 5538 05 


0.0753 258 


7.53 


2 2 P - 


2 2 P 


0.2820 1654 80 


0.0795 333 


7.95 


2 2 P - 


3 2 P 


0.7926 0524 74 


0.6282 231 


62.82 


2 2 P - 


4 2 P 


0.1664 2355 98 


0.0276 968 


2.77 


2 2 P - 


5 2 P 


0.0044 8810 06 


0.0000 201 


0.00 


Total 






0.8107 991 


81.08 


p. 

1 ion 






0.1892009 


18.92 


3 2 P - 


]2p 


0.1602 9051 12 


0.0256 930 


2.57 


3 2 P - 


2 2 P 


0.1587 5209 09 


0.0252 022 


2.52 


3 2 P - 


3 2 P 


0.0646 4956 75 


0.0004 180 


0.04 


3 2 P - 


4 2 P 


0.1782 7236 64 


0.0317 810 


3.18 


3 2 P - 


5 2 P 


0.4245 0682 17 


0.1802 060 


18.02 


Total 






0.2633 002 


26.33 


p. 

1 ion 






0.7366 998 


73.67 


4 2 P - 


1 2 P 


0.1013 4091 61 


0.0102 700 


1.03 


4 2 P - 


2 2 P 


0.1007 7182 29 


0.0101 550 


1.02 


4 2 P - 


3 2 P 


0.1446 2244 15 


0.0209 157 


2.09 


4 2 P - 


4 2 P 


0.0073 1346 56 


0.0000 535 


0.00 


4 2 P - 


5 2 P 


0.1265 8442 00 


0.0160 236 


1.60 


Total 






0.0574 178 


5.74 


p. 

1 ion 






0.9425 822 


94.26 


5 2 P - 


1 2 P 


0.0810 5021 14 


0.0065 691 


0.66 


5 2 P - 


2 2 P 


0.0797 1490 95 


0.0006 354 


0.06 


5 2 P - 


3 2 P 


0.0278 0955 64 


0.0007 734 


2.98 


5 2 P - 


4 2 P 


0.1425 1492 45 


0.0203 105 


2.03 


5 2 P - 


5 2 P 


0.3507 2668 43 


0.1230 092 


12.30 


Total 






0.1605 397 


16.05 


p. 

1 ion 






0.8394 603 


83.94 



°The energies of the incident P-states of Li atom are: 1 2 P -7.4101 3402 79 a.u., 2 2 P -7.3370 5244 47 a.u., 
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3 2 P -7.3117 1131 29 a.u., 4 2 P -7.2996 7937 20 a.u. and 5 2 P -7.2939 6712 18 a.u. For more details of the 
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TABLE V: Transition probabilities between states of D-symmetry for the nuclear (3 -decay of the 
Li atom" to the Be + ion b . 



States Li — > 


Be+ Amplitude 


Probability 


Pi in % 


1 2 D -> 


1 2 D 


0.6142 8174 56 


0.3773 421 


37.73 


1 2 D -> 2 2 D 


0.6684 9584 40 


0.4468 867 


44.69 


1 2 D -> 3 2 D 


0.1470 1259 45 


0.0216 127 


2.16 


1 2 D -> 4 2 D 


0.1666 2545 25 


0.0277 640 


2.77 


1 2 D -> 5 2 D 


0.0801 1513 53 


0.0064 184 


0.64 


Total 






0.8800 239 


88.00 


p. 






0.1199 761 


12.00 


2 2 D -> 


1 2 D 


0.3211 3958 42 


0.3209 154 


32.09 


2 2 D -> 


2 2 D 


0.0738 9620 29 


0.0054 606 


0.55 


2 2 D -> 3 2 D 


0.1751 7333 89 


0.0306 857 


3.07 


2 2 D -> 


4 2 D 


0.0085 5473 79 


0.0000 732 


0.00 


2 2 D -> 


5 2 D 


0.1192 6222 20 


0.0142 235 


1.42 


Total 






0.3713 584 


37.14 


p. 

J ion 






0.6286 416 


62.86 


3 2 D ->• 


1 2 D 


0.4996 1439 98 


0.2496 145 


24.96 


3 2 D -!■ 


2 2 D 


0.5266 3216 24 


0.2773 414 


27.73 


3 2 D -)■ 


3 2 D 


0.1994 4115 14 


0.0397 768 


3.98 


3 2 D -> 


4 2 D 


0.1442 8272 30 


0.0208 175 


2.08 


3 2 D -> 


5 2 D 


0.1192 6222 20 


0.0142 235 


1.42 


Total 






0.6017 737 


60.18 


p. 






0.3982 263 


39.82 


4 2 D -> 


1 2 D 


0.0742 7513 33 


0.0055 168 


0.55 


4 2 D -> 


2 2 D 


0.1217 8600 46 


0.0148 318 


1.48 


4 2 D -> 3 2 D 


0.5611 2979 31 


0.3148 666 


31.49 


4 2 D -> 4 2 D 


0.2051 1260 05 


0.0420 712 


4.21 


4 2 D -> 


5 2 D 


0.0562 9905 16 


0.0031 696 


0.32 


Total 






0.3804 560 


38.05 


p. 

J ion 






0.6195 44 


61.95 


5 2 D -!■ 


1 2 D 


0.0760 3320 23 


0.0058 266 


0.58 


5 2 D -> 


2 2 D 


0.0710 4274 30 


0.0050 471 


0.50 


5 2 D -> 3 2 D 


0.5839 9688 81 


0.3410 524 


34.10 


5 2 D -> 4 2 D 


0.2109 1540 90 


0.0444 853 


4.45 


5 2 D -> 5 2 D 


0.1205 7441 46 


0.0145 382 


1.45 


Total 






0.4109 496 


41.09 


p. 

J ion 






0.5890 504 


58.90 



a The energies of the incident D-statcs of Li atoms are: 1 2 D -7.3355 Of 64 04 a.u, 2 2 D -7.3111 6142 89 a.u., 

24 

3 2 D -7.2331 8155 48 a.u., 4 2 D -6.5887 9667 72 a.u. and 5 2 D -6.6010 1966 62 a.u. 
^The energies of the final states of the Be+ ion arc: f 2 D -13.8780 0469 69 a.u. 2 2 D -13.7793 2881 02 a.u., 



TABLE VI: Study of the accuracy of the wave functions and precision of the calculations for the 
determination of the /3-decay transition probabilities Li atom — > Be + ion. Convergency of the CI 
calculations with respect to the basis. Comparison with Hy-CI calculations. Energy in a.u. 



State 


Basis 


N Amplitude Probability 


Virial 


Energy (CI) 


Ref. Ener." 


Din* DifF in % 


l 2 SLi 


ni 


406 






2.00004 


-7.4748 3857 




3222 




l 2 SLi 


n5 


690 






2.00000 


-7.4763 0444 




1756 




l 2 SLi 


n6 


1388 






2.00000 


-7.4769 2455 




1136 




l 2 SLi 


n7 


2361 






1.99999 


-7.4771 9160 


-7.4780 6032 39 


869 




1 2 S Be+ 


ni 


406 


0.73860 


0.54553 


1.99999 


-14.3190 6528 




5698 


-3.16 


1 2 S Be+ 


n5 


690 


0.77531 


0.60110 


1.99999 


-14.3227 7416 




1989 


2.40 


1 2 S Be+ 


n(i 


1388 


0.72297 


0.52269 


1.99999 


-14.3234 6780 




1295 


-5.44 


1 2 S Be+ 


n7 


2361 


0.77092 


0.59432 


1.99999 


-14.3237 6855 


-14.3247 6317 68 


995 


1.72 


2 2 S Be+ 


n4 


406 


0.51835 


0.26869 


2.00060 


-13.9190 5228 




3737 


0.35 


2 2 S Be+ 


n5 


690 


0.52015 


0.27056 


2.00011 


-13.9207 7083 




2018 


0.55 


2 2 S Be+ 


n(i 


1388 


0.65454 


0.42842 


2.00003 


-13.9215 2875 




1261 


16.33 


2 2 S Be+ 


n7 


2361 


0.63282 


0.40046 


2.00002 


-13.9218 2982 


-13.9227 8926 85 


959 


13.54 


3 2 S Be+ 


ni 


406 


0.08550 


0.00731 


2.00489 


-13.7950 2353 




3693 


0.16 


3 2 S Be+ 


n5 


690 


0.07722 


0.00596 


2.00216 


-13.7966 3281 




2084 


0.02 


3 2 S Be+ 


n(i 


1388 


0.06246 


0.00390 


2.00042 


-13.7974 4352 




1273 


-0.16 


3 2 S Be+ 


n7 


2361 


0.13474 


0.01815 


2.00009 


-13.7977 5382 


-13.7987 1657 


963 


1.25 



"Comparison of energy results are detailed in Ref. [14 1 

6 Difference energy between CI and reference energy in millihartrees. 

c Difference between the transition probabilities calculated by Hy-CI and CI methods in %. 
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